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The first law of thermodynamics states that the average total energy current between different 
reservoirs vanishes at large times. In this note we examine this fact at the level of the full statistics 
of two times measurement protocols also known as the Full Counting Statistics. Under very general 
conditions, we establish a tight form of the first law asserting that the fluctuations of the total 
energy current computed from the energy variation distribution are exponentially suppressed in the 
large time limit. We illustrate this general result using two examples: the Anderson impurity model 
and a 2D spin lattice model. 
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Recent technical advances in the control of nanoscale 
systems have enabled the experimental study of out of 
equilibrium thermodynamics in the quantum regime Hi¬ 
ll]. These new experiments allow for the assessment of 
fluctuations in addition to the mean heat and particle 
currents, thus leading to a renewed theoretical investiga¬ 
tion of the related quantum thermodynamic laws. 

One of the basic questions in this context concerns the 
energy flow between two initially isolated large systems 
A and B. The purpose of this note is to study some con¬ 
sequences of energy conservation on the statistical prop¬ 
erties of this flow. 

By the first law, the average work performed by the in¬ 
teraction coupling the two systems is equal to the average 
heating of the combined system: 

AkUt = AQt = AQf + AQf. 


In the case of a sudden switching on/off of the interaction 
V, the average heating is given by 


AQt = -(V}t + (V}o (1) 


where (• )i denotes the expectation with respect to a suit¬ 
able system state at time t. Whenever V is bounded, Q 
gives 


lim 


AQt 

t 


= 0 . 


( 2 ) 


The individual energy currents AQ^^^ jt are also ex¬ 
pected to reach steady values . They satisfy = 
— J^, and are non-vanishing for systems out of equilib¬ 
rium. 

This note concerns the statistical character of the first 
law related to the thermodynamics of open quantum sys¬ 
tems at the mesoscopic scale. Our main result is a re¬ 
finement of relation ([^. It states that the fluctuations of 
the total energy current are exponentially suppressed in 
the large time limit. 


The nature of work in quantum physics is more sub¬ 
tle than in classical physics [10]. In the 1990’s Lesovik 
and Levitov introduced the concept of the Full Counting 
Statistics (FCS) in the study of charge transport [TT] . 
The use of the FCS in the definition of work in quantum 
physics appeared in the early 2000’s in the works of J. 
Kurchan and H. Tasaki on the extension of the fluctua¬ 
tion relations to quantum systems [nun]. The emerg¬ 
ing idea is that in quantum mechanics work should not 
be understood as an observable. Instead, the work per¬ 
formed during a given time period is identified with the 
energy variation AE observed in a repeated measure¬ 
ment protocol where the system energy is measured at 
the beginning and at the end of the period. The distri¬ 
bution of the measured energy variation, Ft{AE), is the 
work FCS (we comment on terminology in footnote [14] 1. 
This change of perspective opened a whole new area of 
research da m]. In particular, it allowed for the ex¬ 
tension of the fluctuation relations to quantum systems 

dllllllMH]- 

The fluctuation relations are intimately related to the 
second law of thermodynamics and have been extensively 
studied dll El . Regarding the first law, the well 

known identity 


Et{AE) = AQt 


and ([^ give 


(^) = “ W 

where Et denotes the expectation with respect to the 
FCS distribution dHl E]- In this note we sharpen 
([^ by showing that, under very general conditions, the 
exponential moment 
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remains bounded as t —>■ oo where the constant am > 0 
is a measure of the regularity of the interaction V (see 
([^ below). 

Until recently, the first law and energy conservation in 
the FCS setting have received little attention in the liter¬ 
ature. In the case where A and B are thermal reservoirs, 
the FCS of the total energy current was previously stud¬ 
ied theoretically in m- The works [231 concern the 
FCS of energy transfer in the thermalization process of 
a finite level quantum system in contact with a thermal 
bath, a problem which is radically different from the one 
considered here. We also emphasize that here we are only 
interested in the FCS of the total energy, and not in the 
FCS of the individual energy variations . 

We start with a system described by a finite dimen¬ 
sional Hilbert space where the superscript L refers 
to the size of the system. Taking L —>■ oo corresponds 
to the thermodynamic limit. The limiting objects will be 
denoted without the superscript. Let 
be the Hamiltonian of the joint but non-interacting sys¬ 
tem A+B. The evolution between the two measurements 
of is generated by Hy^ = where 

denotes the interaction coupling A and B. The initial 
state is described by the density matrix 

Let denote the projection on the eigenspace asso¬ 
ciated to the eigenvalue e in the spectrum The 

measurement of at initial time t — 0 gives e with 
probability tr(Pi^^p*^^)). After the measurement the sys¬ 
tem is in the projected state 

pWpWpiL) p{L)y 

The second measurement of at a later time t gives 
e' with probability 

tr /tr(pWp(^)). 

It follows that the probability of observing the energy 
variation AE in this measurement protocol is 

f[^\AE)= tr 

e' — e—AE 

The moment generating function of the Full Counting 
Statistics is 

x['^\a)= [ e“^^dp(^)(AP) 

Jr 

=tr e*‘^-’) 

where 

= E 

We assume that for a purely imaginary, the limit 

\im = Xtia) (4) 

L—¥oo 


exists and is a continuous function of a. This assumption 
is harmless and easy to verify in most concrete models of 
physical interest. By Levy’s continuity theorem m, @ 
implies that the thermodynamic limit lim^^-^oo = Pt 
exists. The probability distribution is the FCS of the 
thermodynamic system. 

Let 

=2\a\ max || 

-1<S<1 

and 

R{a) = sup R^^\a). 

L 

Note that R{a) takes values in [0, oo] and is an even func¬ 
tion. Moreover, R[a) > R{a') if a > a' > 0. Our regu¬ 
larity condition is that there exists am > 0 such that 

R{ ) < 00. (5) 

We emphasize that ([^ is the only regularity assumption 
we require and that no further hypothesis on the dynam¬ 
ical behaviour of the system is needed. We also make no 
assumptions on the initial state of the system. 

Our main result is the following strengthening of ([^: 
Theorem For alH > 0, 

Et (e“™|A^I^ < (6) 

An immediate consequence of this result and Cheby- 

shev’s inequality [ 23 ] is that for any e > 0, 

Pt (7) 

Note that if R{a) < oo for all a, then 

2^R{CP)-Ct (8) 

for any C > 0. 

The estimates Q and ^ can be interpreted in terms 
of the large deviation theory |2S| (see [23]). For example, 
Q implies that the large deviation rate function of the 
random variable |AP|/t satisfies I{s) = oo for s 7 ^ 0, 
and that the large deviations are completely suppressed 
in the large time limit. 

The main novelty of our proof is the derivation of a 
time independent bound for Xt^'^ inspired by the bounds 
proposed in [2T]. The derivation is based on two well- 
known inequalities. The first is 

tr(Ay) < ||A||tr(y) 

which holds for any two non-negative matrices X, Y. The 
second states that for any two self-adjoint matrices T, S, 

||gT-|-Se-T|| ^ gmaxo<,<i ||e=’^Se“=’^|| ^ 


(9) 
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To prove this inequality, let r(s) = e^^-^+'^^e Then 
one has 

5,r(s) = r(s)e"^5e""^, r(o) = I 

Using 

\\dsT{s)\\<\\ris)\\\\e^^Se-^^\\ 

and Gronwall’s inequality we obtain The bound ([^ 
is similar but unrelated to the bound ^.10) of [27]. 

The proof of Q proceeds as follows. For a real we set 


theorem (see Appendix B in [TS]) give that for all com¬ 
plex a with real part Re(a) in the limit 

limL_>oo = Xtio:) exists. Moreover, for such a’s, 

Xtia) = [ e“^^dPt(AU) 

Jr 

and 

It follows that 


A = e 




(L) 

V e 




(L) 


|Xt(±am)| < 


and 


The last estimate gives 


y — 0 0 2 0 2 Q 2 0 2 0^‘'‘'‘^ 

(note that and commute). Observe that 


{L) 

V 


xi^Ha)=tr(XY) 

and that X, Y are non-negative matrices. We then use 
the first inequality to derive the estimate 

x[^\a) < ||A|jtr(r) 


where 

||A|| = 

and 

tr(r) = tr (et<’e-t"'"’pWe-t^<">et<’) . 

The cyclicity of the trace gives 

tr(r) = tr (e-t"‘"’e“<’e-t^'"’pW) . 

Applying the first inequality once again and using that 
tr(p(^)) = 1, we derive 


tr(F) < ||e 2 -'^v 


(L) 




(i) , 


r(i) 


Hence 

Using the second inequality with 


0 2 V 0 2 ^ 


E( < Xti-am) + Xt{am) < 

and the theorem follows. 

a. Spin-fermion models. Electronic transport 
through a ID-lattice containing a single magnetic impu¬ 
rity is a typical problem involving bounded interactions. 
The Anderson model [2^ [29| commonly used to study 
this question is a specific example of a general class of 
spin-fermion models to which our main theorem applies. 

The study of the FCS of charge transport through the 
impurities in such models is an active field of research [SUF 
134] . We emphasize, however, that here we are only con¬ 
cerned with the statistics of the total energy. 

The impurity is described by a quantum dot support¬ 
ing four different eigenstates: empty, occupied by a single 
electron with either spin up or spin down, or occupied by 
two electrons with opposite spins. The remaining parts 
of the lattice, regarded as fermionic (say left and right) 
reservoirs at different chemical potentials, are described 
in the tight binding approximation. 

Here, the subsystem A is the left side of the lattice 
together with the impurity. The lattice right side is the 
subsystem B. 

The operator c*^^^{x) {ci/r^cr{x)) creates (annihilates) 
an electron with spin a at the lattice site x of the left (x < 
0)/right (x > 0) reservoir. Similarly, the operator d* (da) 
creates (annihilates) and electron with spin a in the dot. 
The anti-commutation relations ^,(x'y\ = 

dx,x'da,a' and {dcr,ci*/} = 6a,a' hold while the c operators 
commute with the d operators. We use the shorthand 
ci/r,(y{4') = The reservoh Hamiltonians 

are 


we obtain 






The regularity assumption ([^, the existence of the limit 
Q for purely imaginary a’s, and Vitali’s convergence 


Hi= cl^{x)cp^{x'), 

cr=:±;x,x'<0 
|x-x'| = l 

with a similar expression for H^.. Let be the discrete 
Laplacian of the left/right part of the lattice. Since 
is a bounded operator, 
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for all real a. In particular, for all a, 

< oo. (10) 

The total Hamiltonian is 


H = Hs +Hi + Hr 

where Hs = e d’^d^+Ud*_^d+d‘^d- is the Hamiltonian 
of the dot. Regarding the subdivision in A/B subsystem, 
we have Ha = Hi + Hg and Hb = Hr- The coupling of 
the conduction electrons with the dot is described by 

R — ^ ^ {,dfy{ci^(y{yi^cj^ T Cr^a “b h.C.) 


for some coupling functions vifr,a{x)- In the context of 
the Anderson model, the superscript L refers to the con¬ 
finement of the reservoirs to the finite part of the lattice 
defined by |a;| < L. Such confinement is necessary to 
allow for a meaningful definition of the repeated mea¬ 
surement protocol leading to the FCS. The limit L —^ oo 
restores the extended reservoirs. It follows from relation 


(101 that R{a) is finite for all a, and that our theorem 
holds for all am > 0. Hence we have inequality ®: 




for any e > 0 and any C > 0. 

We also note that one can consider instead the FCS of 
H' = Hi + Hr by setting V = Hs + V. Then H^ = Hi 
and H'b = Hr- One then obtains the same result by 
replacing AA with AA'. The energy of the impurity is 
irrelevant in the large time limit. 

b. Spin systems. Another popular class of models 
involving bounded interactions are locally interacting 
spin systems. In |26j we prove that, under general con¬ 
ditions, our theorem applies to locally interacting spin 
systems in arbitrary dimension. Moreover, for ID spin 
systems with finite range interactions, Araki’s results |85j 
give that R{a) < oo for all a, and hence that our the¬ 
orem holds for all am > 0. We restrict ourselves to the 
description of a simple example. 

Consider a 2D square lattice of ^-spins. Let A^ C 
be the finite sub-lattice of size 2L x 2L. We denote by 
A^ its left/right half. Subsystems A and B are the spins 
in A£ and Aj respectively (see Figure]^. 

The system Hilbert space is RH) = 

Hamiltonian is that of an XY-spin model where the spins 
on A/ do not interact with that on A/( [HH] : 

hH) ^ 

with 

x,y nearest 
neighbors in a/ 


A B 



FIG. 1. A partitioned finite spin system A A- B. Solid lines 
represent the nearest neighbour coupling J and dashed lines 
the interaction Kr,y between the 2 subsystems. 


where J is a coupling constant. The interaction is 

xeA'[,yeA+ 

where 


if a; = ( 0 ,^ 2 ) G A/ and y = ( 1 ,X 2 ) S Aj and = 0 
otherwise. The boundary between the two halves of the 
lattice is between the lines xi = 0 and xi = 1. Note that 
the interaction intensity decreases as one moves away 
from (0,0). An assumption of this type is necessary if 
is to remain bounded in the thermodynamic limit 
L —y 00 . 

For this model one can show that there exists am > 0 
such that (§ holds and that our theorem applies. Hence 
we have inequality Q: 

Ft - 2e"‘“’"+^(“™) 

for any e > 0. 

c. Discussion. Under a general condition on the reg¬ 
ularity of the interaction evolution in imaginary time, we 
have proven a sharp form of the first law of thermody¬ 
namics for the FCS of energy variation. 

Our result holds for any initial state of the system. 
If one assumes that systems A and B are initially in 
thermal equilibrium at temperatures Ta and Tb, then 
the suppression of the fluctuations of the total energy 
current can be also proven by following the arguments of 

Under additional assumptions it is possible to deal with 
cases where several reservoirs drive the joint system to¬ 
wards a non-equilibrium steady state and to derive prop¬ 
erties of the joint distribution of the energy variations 
in each part of the system. A more strict condition on 
R{a) allows for the generalization of a symmetry of the 
limiting cumulant generating function proposed in |21j . 
Combined with time reversal invariance this leads to On- 
sager’s reciprocity relations. We investigate these topics 
in 
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In the present note we have limited ourselves to 
bounded interactions. The case of unbounded interac¬ 
tions (an example is the spin-boson model) is more tech¬ 
nical and requires a separate analysis based on an appli¬ 
cation of Ruelle’s quantum transfer operators [H]. Al¬ 
though the physical picture emerging from this analysis is 
of an independent interest, the final results are much less 
general than in the case of bounded interactions mm- 
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